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ABSTRACT 

Previous observations show that the hypercompact HII region W51e2 is sur- 
rounded by a massive molecular accretion flow centered on the HII region. New 
observations of the H53a radio recombination line made with the VLA at 0.45 
arc second angular resolution show a velocity gradient in the ionized gas within 
the HII region of > 500 kms— 1 pc _1 comparable to the velocity gradient seen 
in the molecular accretion flow. New CO line observations made with the SMA 
at arc second angular resolution detect a molecular bipolar outflow immediately 
around the W51e2 HII region and extending along the axis of rotation of the 
molecular flow. These observations are consistent with an evolutionary phase for 
high mass star formation in which a newly formed massive star first begins to 
ionize its surroundings including its own accretion flow. 

Subject headings: HII Regions 



Introduction 



1.1. HII regions within star-forming accretion flows 



Molecular line observations of several massive-s tar forming regions show evidence for 



star-forming accretion flo ws around small HII regions (IGuilloteau et al 
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G10.6-0.4, W51e2, G28.20-0.04, and G24.78+0.08) are of hypercompact (HC) or ultracom- 
pact (UC) size and are bright enough to require an ionizing flux at a level produced only 
by one or more O-type stars. Thus the presence and brightness of the HII regions at the 
centers of these flows indicate that the accretion is associated with the formation of the most 
massive stars. 

The presence of the HII regions within molecular accretion flows poses some simple yet 
interesting questions. How does the development of an HII region affect accretion and the 
evolution of a young star? Does the formation of an HII region prevent further accretion 
and growth of the star? How does an HII region affect a pre-existing, bipolar outflow (BPO) 
of the type that is typically associated with the star-forming accretion flows of lower mass 
stars? 

The interaction of HII regions with accretion flows and associated BPOs occurs in an 
evolutionary phase in massive star formation common to all young O stars. Stellar structure 
calculations indicate that massive stars with accreti on rates simil ar to those observed around 
the HII regions listed above (1(T 3 - KT 2 M yr _1 , dKetoll2002al table 1)) should b erii 
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nuclear burning at a mass well below that of an O star (Stahler. Shu fc Taam 198ol. 1981a bj; 



Beech fc Mitalas 111994 : iNorberg fc Maeder 112000 : iBehrend fc Maeder 



20061 ). These calculations also indicate that the temperature and lumi- 



nosity of the accreting star is about the same as a ZAMS star of equivalent mass. Thus an 
accreting star begins to produce an HII region once the star reaches a mass, temperature and 
luminosity equivalent to early B, about 15 - 20 M Q . If the star is to accrete the additional 
mass to become an O star, then the accretion must continue past the mass when an HII 
region to forms. Thus the interaction of accretion flows and HII regions is inevitable in the 
early evolution of massive stars; however, this interaction is currently not well understood. 

Since bipolar outflows are invariably associated with star-forming accretion flows of lower 
mass stars up to later-type B, an HII region that develops within a pre-existing molecular 
accretion flow must also interact with the molecular bipolar outflow that we expect to be 
associated with the accretion flow. There are few detections of outflows from young O stars. 
Young O stars are rare both because of the steep stellar mass spectrum and because of the 
short time that they spend in the accretion phase. Furthermore, they are difficult to identify 
because they are deeply embedded in molecular clouds and always at kpc distances. Maybe 
we have not yet found enough examples, or maybe the HII regions that surround all young 
O stars modify or terminate the outflows. 
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1.2. The ionized accretion flows 



Of the HII regions associated with molecular accretion flows, three have been observed 
in radio recombination lines (RRL) with high enough angular resolution to map the velocities 
of the ionized gas within the HII regions. All three HII regions show rotation in the ionized 
gas and one also shows infall. The similarities of the velocities in the molecular and ionized 
phases suggest that the flow is continuous across the ionization boundary. The UCHII region 
G10.6-0.04 which contains a small group of O stars is the best studied example, and t he one 



that also shows an inward flow of the ionized gas fcetoll2002al : iKeto fc Woodl hood) . The 



G34. 3+0.2 HII region shows rotation in its cometary shaped ionized outflow (IGarav et al. 



1986h that is perpendic ular to a rota t ionall y flattened molecular accretion flow (IKeto et al.l 
1987bl ). More recently, ISewilo et al.l (120071 ) report evidence for rotation of the ionized gas 



within the G28.20-0.04 HII region in the same orientation and velocities as the rotation 
seen in the "97 kms" 1 " component of the NHs(3,3) line previously mapped by ISollins et al. 



( l2005bl ). In addition to t hese three, rotation of the ionized gas within the G45.07+0.13 HII 
region is also reported in IGarav et al.l (119861 ) , but observations of the s urrounding molecular 



gas are inconclusiv e as to the presence of a molecular accretion flow (jCesaroni et al.l 11994 
Wilner et al1ll996h . 



Inward flowing ionized gas has been observed only in the large, cluster-scale HII region 
G10.6-0.4. This HII region is the only one of the three observed examples that has sufficient 
stellar mass inside the HII region (a few hundred M ) that the Bondi-Parker transonic radius 
R g = GM/2c 2 s is large enough (several thousand AU) to be observationally resolved at at the 
distances (6-8 kpc) of these three HII regions. The other two HII regions might also have 
inward flowing ionized gas, but if they contain only one O star or even a binary or triple, 
then the inward flow will be on a length scale that is too small to resolve observationally. 



1.3. The massive bipolar outflows 



There have been a few single di sh telescope surveys tha t stud y B POs in star forming 
region s with UC and HC HII regions. IShepherd fc Churchwelll (119961 ) and lKlaassen fc Wilson 
( 2007 ) looked for BPOs in massive star-forming regions with UC/HC HII regions. IShirley et al 
( 120031 ) selected star-forming regions with water maser emission. Because water masers are 
associated with high mass stars, these regions also often contain UC or HC HII regions. 
These surveys detect BPOs with mass outflow rates and energies that indicate that they are 
driven by massive stars. However, these single dish observations do not have the angular res- 
olution necessary to determine whether the BPOs are associated with the bright HII regions 
and therefore O stars or with other nearby massive stars currently equivalent to ZAMS B 
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without HII regions. In order to determine whether a BPO is associated with a particular 
accreting massive star or an HII region, we need higher resolution spectral line observations 
around the HII region in molecules that best trace bipolar outflows. 



New observations 



We chos e to observe the W51e2 HII region because of its well-studied molecular ac- 



cretion flow 



Sollins et al. 



Rudo 



2004 



ph et al 



1990: iHo fc Youndll996l : IZhang fc Holll997l : IZhang et al.lll998 



Young et al.lll998l ). In this paper we report on VLA observations of the 



H53a radio recombination line (RRL) and on SMA observations of CO(J=2-l). We observed 
the highest frequency RRL accessible by the VLA because the higher frequency lines have 
lower optical depth, and we can observe the velocities of the gas deeper inside the HII region. 
We observed the CO in a transition above the ground state and at ~ 1" angular resolution 
to improve the likelihood of observing a BPO associated with hot molecular gas near the 
HII region. 



2.1. Experimental setup 

Our H53a observations of the ionized gas in W51e2 were made at the National Radio 
Astronomy Observatory Very Large Array (VLA)0. The H53a line was observed in the C 
array at a resolution of 0.45". The 12 CO(J=2-l) line was observed at the Submillimeter 
Array (SMA) @ in its extended configuration, resulting in a spatial resolution of ~ 1". The 
C O line was detected simultane ously with H30a, and both sets of observations are described 
Keto. Zhang fc Kurtz ( 2008 ). The James Clerk Maxwell Telescop^fl was used to observe 



in 



the larger scale structure of the molecular gas, and these observations were combined with 
the SMA observations using the non-linear image combination program MOSMEM in the 
MIRIAD data reduction package (as described in Klaassen et al. in prep). 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

2 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 

3 The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 
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2.2. The ionized flow 



Figure 1 shows the average velocity of the H53a line. The velocity gradient is approxi- 
mately 8 kms -1 between 57 and 65 kms -1 VLSR with a position angle —30° (west of North). 
This velocity pattern is consistent with either rotation or outflow. We prefer the interpreta- 
tion that the H53a velocities indicate rotation for a couple of reasons. First, the observed 
velocity gradient is perpendicular to the CO outflow. Bipolar outflows associated with lower 
mass stars are always oriented approximately perpendicular to the rotation. Second, the 
interpretation of rotation is consistent with the magnitude and d i rection of the rota tional 
velocity gradient observed in NH 3 and CH 3 CN (IZhang fc Holll997l ; IZhang et al.lll998l ). 



The clearest signature of rotation in the molecular gas is in figure 7 of IZhang fc Ho 



(119971 ) in the position- velocity diagram of NH3(3,3) along a line at position angle of 135°, 
and consistent (135 - 180 = —55°) with the position-angle of the velocity gradient in the 
ionized gas (—30°). We prefer to use the NH 3 (3,3) line in this interpretation rather than 
the (2,2) line because the higher excitation line should derive from gas closer to the hot 
HII region. However, there is considerable uncertainty in locating the exact angle of the 
axis of rotation using the molecular line observations. These observations map the flow at 
relatively large scales before the accretion has spun-up by angular momentum conservation 
to higher rotational velocities. The observed infall velocities are comparable or greater 
than the rotational velocities, and the rotational component is difficult to extract. For 
example, comparison of the position angle (-55°) of the rotational gradie nt derived from the 
NH3 observations with the position angle (—70°) derived from CH3CN ( IZhang et al.lll998l ) 
suggests an uncertainty of several tens of degrees. 

Our observations do not fu l ly res olve the HII region, but the 1.3 cm continuum obser- 



estimated by Zhang fc Hoi ( 1997 ) and 



vation of lKeto. Zhang fc Kurtzl (12008) indi c ates a FWHM of 0.4" similar to the upper limit 



Scottl (119781 ). The radius, R, associated with the veloc- 



ity difference AV, is uncertain because of our low angular resolution. So we use the FWHM 
of the continuum emission as a characteristic size. Assuming that the observed velocities are 
rotational, the mass of the star within the W51e2 region is estimated as M = AV 2 R/2G. 
Assuming a distance of 8 kpc, we derive a lower limit to the velocity gradient of > 500 kms -1 
pc _1 and a dynamical mass of V 2 R/2G > 15 M This mass is roughly consistent with a 
previous est imate of about 20 M Q derived from the emission measure of the radio continuum 
JScottlll978h . 



This estimated ve l ocity gradient in the ionized gas is similar to that measured in molec- 



ular gas. IZhang fc Hoi (119971 ) measure a velocity gradient in NH3 of 500 kms 1 pc 1 centered 



around 58 kms 1 VLSR. Owing to limi ted angular r e soluti on, the exact value of the velocity 
gradient may be somewhat uncertain. IZhang et al.l (119981 ) derive a lower value (100 kms -1 
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pc _1 ) for the velocity gradient from observations of CH3CN. Nonetheless, the similarity in 
the magnitude and direction of the velocity gradients observed in the molecular and ionized 
gas suggests that there is a continuous flow of gas from the molecular to the ionized phase 
such that the ionized gas retains some memory of the velocities of the molecular accretion 
flow. 



2.3. The bipolar outflow 

Figure 2 shows the bipolar outflow in CO. The orientation of the flow is along the 
rotation axis of the ionized and mo lecular accretion flows defined by the H53a RRL and NH 3 



observations of (jZhang fc Holll997l ). The outflow is seen as blue and red lobes in the maps of 
integrated emission. In figure 2, the blue lobe is shown in the emission integrated from 20 to 
44 kms -1 (blue) and the red lobe from 68 to 92 kms -1 (red). These ranges are both 24 kms -1 
wide, starting at velocities ±12 km s -1 from the VLSR of the larger scale molecular gas (56 
km s" 1 ). The blue shifted lobe is much brighter than the red one and better collimated. 
The better collimation of the blue flow results in a smaller volume. If the entrained mass is 
estimated from the volume and brightness, the two flows have approximately the same mass 
(3-4 M ) and kinetic energy (~ 3 x 10 43 J), as calculated from zeroth and first moment 
maps of the outflow (Klaassen et al. in prep). 



3. Implications 



A conceptual model for the ionized flow in W51e2 may be constructed as a composite of 
a couple of model calculations in the literature. The expected structure (morphology) of an 
HII region that dev elops withi n a large scale molecular accretion flow may be calculated from 



ionization balance (lKetoll20071 ). In particular, figure lc of that paper shows a nearly spherical 
HII region around a rotationally flattened molecular flow that may be relevant for W51e2. 
These calculations describe the structures of HII regions but do not calculate the velocities 
or dens ities of the io nized gas. Rather they assume an ionized flow based on the Parker wind 



model (jParkerlll958l ). The velocities and densities of the gas that flows into an HII from a 



disks (PEDs) ( 


Hollenbach et al. 


1994; 


Yorke 


1995; 


Yorke & Welz 


1996: 


Lizano et al. 


1996; 


Johnstone. Hollenbach. & Ballv 


1998: 


Lueo. Lizano & Garav 


2004 


). However, these models 



do not calculate the structure of the molecular accretion flow and HII region for the types 
of molecular accretion flows that are observed around HII regions. Rather they assume the 
accretion flow is in the form of an infinitely thin disk. In the absence of a more complete 
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calculation we may conceptually combine the flow morphology from iKetol (120071 ) with the 
velocities and densities from the PED models to interpret the the observations of W51e2. 
The spin-up of the flow seen in the molecular line observations suggests that the large scale 
molecular accretion flow becomes progressively more flattened at smaller radii. If the flow 
is sufficiently dense at its mid-pla ne, it is able to resist ionization and extend into the HII 
region as described in IKetol (120071 ). Across the top and bottom boundaries of the flattened 
molecular flow there is an ionization front that continuously supplies gas to the HII region as 
described in the PED models. Since there is no obvious force to immediately stop the rotation 
of the gas as it crosses the ionization front, the ionized gas is also rotating. Observations 
of high frequency RRL such as the H53a are dominated by the densest gas that is near the 
ionization front. The velocities of the H53a line thus show the rotation in the ionized gas 
coming off the rotating molecular accretion flow. 

The density of the ionized gas decreases rapidly as it accelerates off the ionization 
front around the flattened molecular accretion flow. The exact rate of the density decrease 
depends on the geometry of the flow. The Parker wind model assumes a purely spherical flow 
while the PED models describe the flow off a flat, rotating disk. In both cases, the density 
distribution is approximately exponential (hydrostatic) where the flow velocities are subsonic 
and closer to a power law past the transonic point. Observations in the radio continuum 
should thus measure a rising spectral energy distribution (SED) consistent with the steep 
density gradient in the outflowing ionized wind. In the case of W51e2, the observed SED 
may be reproduced by a mod el with a density gradient scaling with radius as the —2.5 power 
fceto. Zhang fc Kurtz! [io08h . 



Th e PED mode ls show that the ionized gas flows along the rotation axis, and the struc- 



tures in IKetol (120071 ) show that this flow can be further channeled by the density structure of 
the molecular accretion flow. This outflowing ionized gas could drive a molecular BPO, but 
this BPO would be quite different in collimation and velocity structure from those BPOs 
associated with the accretion flows of lower mass stars. At the moment, the driver of the 
BPO that we see in CO in not known. It may be the original BP O, but now flowing; th rough 
the HII region, or it may be a remnant of a now terminated flow (IKlaassen et al.ll2006l ) or it 
might be driven by the ionized outflow. 



4. Conclusions 

Our new radio recombination line observations of W51e2 detect rotation in the ionized 
gas consistent with an accretion flow, and our new CO line observations detect a molecular 
bipolar outflow along the rotation axis. 
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The observations are consistent with the structure expected at the time that an accreting 
massive star begins to ionize its surroundings including the accretion flow responsible for the 
formation of the star. 
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Fig. 1. — Velocity of the H53a line in W51e2 in color. The color scale ranges from 57 to 65 
kms -1 . The contours show the 7 mm continuum emission at 2, 4, 10, 30, 50, 70 and 90% of 
the peak emission of 0.15 Jy/beam. Coordinates are in the B1950 epoch. 
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W51e2 C0(2-1) Integroted Intensity 




Fig. 2. — Integrated emission of the CO(2-l) line in color. The emission from 20 to 44 
kms -1 is shown in blue and from 68 to 92 kms -1 in red. The blue scale ranges from -17 to 
47 K kms -1 . The red scale ranges from -9 to 16 K kms -1 . The contours show the 1 mm 
continuum emission at 30, 50, 70 and 90% of the peak emission of 2.1 Jy/beam. This map 
of W51e2 is on a larger scale than figure 1. Coordinates are in the B1950 epoch. 



